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SMALL-ANGlf STABILITY ANALYSIS OF A LINEAR CONTROL SYSTEM 
FOR A HIGH- POWER COMMUNICATION SATELLITE 
by Thomas A. O'Malley 
Lewis Research Center 

SUMMARY 

This report presents a small-angle stability analysis for one particular configura- 
tion of a high-power communication satellite having a linear control system. Both the 
central body and the solar array are treated as rigid bodies. The configuration studied 
is one for which the solar array has one degree of freedom with respect to the central 
body, and the antennas are rigidly attached to the central body. The control system 
studied consists of three-axis control of the central body and one-axis control of the 
solar-array rotation relative to the central body. The results of the report yield pre- 
liminary indications of the relation of stability to satellite inertias and control gains. 

Assumptions are made which allow the equations of motion to be linearized and 
which make the resulting stability analysis tractable. The eighth-order characteristic 
equation is separated into two second-order equations for the roll and yaw axes of the 
spacecraft and a fourth-order equation for the pitch axis. For positive control gains, 
stability is achieved for the roll and yaw axes. Not all positive values of control gains 
will yield stability for the pitch axis. Pitch-axis stability is achieved only when the con- 
trol gains are sufficiently large. 

When the central body and solar array do not have coincident centers of mass, 
steady- state pointing errors result. For reasonable values of inertias and control 
gains, the steady-state pointing errors are negligible. The degree of stability is not 
significantly changed when the mass centers are not coincident. 


INTRODUCTION 

The next generation of high-power communication satellites will use large sun- 
oriented solar arrays as a primary power source. Typical configurations of such com- 
munication satellites have one or more antennas mounted on a relatively dense central 



body. The central body must be tightly controlled in order to transmit narrow radio- 
frequency beams from synchronous equatorial orbit to high-gain receiver antennas on 
the ground. The dual control requirements of earth -pointing the central body and sun- 
pointing the solar array require one or perhaps two degrees of freedom between the cen- 
tral body and the solar array. 

To analyze the attitude control of such communication satellites, the equations of 
motion of multibody spacecraft have to be derived. General methods of deriving these 
nonlinear equations have appeared in a number of reports (refs. 1 to 3) in recent years. 
The process of determining the stability of a control system often involves using stand- 
ard techniques, such as root locus, which can be applied only to a system of linear 
equations. This report presents the assumptions and analysis required to reduce the 
nonlinear dynamic equations to a form in which the standard linear stability techniques 
can be easily applied. 

The dynamic equations used in this report are for a spacecraft consisting of two 
rigid bodies - a central body and a solar array. The configuration studied is one for 
which the solar array has one degree of freedom with respect to the central body, and 
the spacecraft antennas are rigidly attached to the central body. The control system 
studied consists of three-axis control of the central body and one-axis control of the 
solar -array rotation relative to the central body. 

The method of attack is to derive the generalized equations of motion. The full 
derivation is given in appendix A. The generalized equations are then linearized in ap- 
pendix B. Symbols are defined in appendix C. The characteristic equation of the sys- 
tem is obtained from the linearized equations of motion. Approximations are made to 
separate the eighth-order characteristic equation into two second-order equations for 
the roll and yaw axes of the spacecraft and a fourth-order equation for the pitch axis. 
Stability is analyzed first for the case of the central body and solar array having coin- 
cident centers of mass. The analysis is then generalized for distinct centers of mass. 


COORDINATE SYSTEMS 

Before deriving equations of motion, it is necessary to define coordinate systems 
suitable to the satellite configuration. The configuration to be studied in this report, 
shown in figure 1, is the one proposed for the Canadian Technology Satellite. The same 
Configuration is also being considered for the ATS-H and ATS-I satellites. 

These satellites have a dense central body with rigidly attached antennas and a large 
solar array. The attitude control system maintains the central body in an earth-oriented 
position in order to satisfy antenna -pointing requirements. The large solar array is 
sun-oriented. To maintain earth-pointing and sun-pointing requirements, the solar ar- 
ray rotates with respect to the central body at a nominal angular velocity of one revolu- 
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tion per day. Since there is only one degree of rotational freedom of the solar array 
with respect to the central body, the nominal sun-pointing error is zero only during the 
vernal and autumnal equinoxes, when the sun is in the equatorial plane. The nominal 
sun-pointing error is a maximum (23. 5°) during the solstices. 

The central-body coordinate system and the orbital reference coordinate system 
are shown in figure 2. The orbital reference coordinate system, specified by unit vec- 
tors x , y , z r along the coordinate axes, has its origin at the center of mass of the 
central body. The positive z r -axis points to the center of the earth, the positive x r ~ 
axis lies along the velocity vector of the assumed circular orbit, and the y r ~axis com- 
pletes the right-handed set. The central-body coordinate system, specified by unit vec- 
tors Xp yp Zj along the coordinate axes, is fixed in the central body and has its origin 
at the central-body center of mass. The x^, y^, and z^ axes are oriented such that 
they coincide with x f , y r , and z r when the central body is in its nominal orientation. 
The Xj, yp and ~z^ axes will be called the roll, pitch, and yaw axes of the central 
body. The central-body coordinate system is obtained from the orbital reference co- 
ordinate system by three rotations: (1) rotation about x r ~axis through an angle 0p 
(2) rotation about the y r ~axis through an angle 0 2 ; (3) rotation about the z r -axis 
through an angle 0 g. The objective of the central-body control system is to null out 
the pointing errors 0j, 0 2 , 0g. 

The solar-array coordinate system, shown in figure 3, is fixed in the solar array 
and is specified by unit vectors x 2 , y 2 , and z 2 along the coordinate axes. The origin 
of the solar-array coordinate system is at the solar-array center of mass. The posi- 
tive y 2 ~axis is along the solar -array shaft and is parallel to the positive y^-axis. The 
z 2 ~axis is perpendicular to the solar -array surface, and the x 2 -axis completes the 
right-handed set. The x 2 , y 2 , and z 2 axes will be called the roll, pitch, and yaw axes 
of the solar array. Let x'p y'j, and z'j be unit vectors along the axes of a coordinate 
system whose axes are parallel to the central-body axes but whose origin is at the solar- 
array center of mass. The angle cp shown in figure 3 is the angle between x'^ and x 2 , 
or equivalently the angle between z'j and z 2 . 

Let 6 be a unit vector originating from the solar-array center of mass and pointing 
to the sun, as shown in figure 4. Then the true solar-array pointing error is the angle 
between z 2 and 6. This pointing error, however, may have a nominal value as large 
as 23. 5°. For attitude control purposes, we assume that the error sensed by the con- 
trol system is the angle o shown in figure 4. The vector 6p is a unit vector obtained 
from the projection of 6 into the x 2 - z 2 plane. The angle subtended by z 2 and 6p 
is CT. 

The four rotational degrees of freedom of the system are specified by the angles 
0p 0 2 , 0g, and cp. The angles 0p 0 2 , and 0g will be small angles; whereas, the 
angle cp nominally increases at the rate of 2n .radians per day. Since the four angular 
coordinates used for the stability analysis must remain small in order to justify linear - 
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izing the equations of motion, cp cannot be used as one of the coordinates. The fourth 
coordinate must be defined as a perturbation from a nominal motion. 

The geometry for defining a suitable fourth coordinate is shown in figure 5. The 
plane of the figure is the equatorial plane, and the line L is the projection of the earth- 
sun line into the equatorial plane. The line L may be considered fixed for the time 
span considered. Let x' r (t) and z^,(t) be unit vectors, at time t, along the axes of a 
coordinate system whose axes are parallel to the orbital reference axes but whose origin 
is at the center of the earth. The angle subtended by the line L and the vector z^.(0) 
is a. The angle subtended by x^.(0) and x^(t) is a) Q t, where cu Q is the angular veloc- 
ity of the synchronous orbit. 

The fourth coordinate 9 4 is defined as 

0 4 (t) = <o(t) - (a + w 0 t) (1) 

If the spacecraft is in its nominal position at time t, the line L is coincident with z 2 (t), 
and z r (t) is coincident with z^t). Thus, cp( t), the angle between z 2 (t) and Zj(t), will 
be the angle between L and z^(t), which by figure 5 is the angle a + a> o t. From equa- 
tion (1), we conclude that 0^(t) is zero when the spacecraft is in its nominal position. 
Thus, 9 ^ is a small angle and is a suitable fourth coordinate. 

When the sun is in the equatorial plane, the solar-array error a is the sum of 0 2 
and 0 4 . When the sun is not in the equatorial plane, the central-body roll and yaw 
errors 9 4 and 9 2 will also contribute to c r. In the general case, a is given by (see 
appendix B) 


a = (G sin + 0 2 + (G cos cp)9^ + 9 4 


( 2 ) 


where G is defined as 


-sin sin i 

G = — 8 (3) 

yi - sin^0 g sin^i 

The angle i of equation (3) is the inclination of the equatorial plane to the ecliptic plane, 
and the angle 9 is the apparent angular displacement of the sun from the autumnal 
equinox measured in the ecliptic plane, as shown in figure 6. 

ASSUMPTIONS 

It is necessary to make assumptions which will justify linearizing the equations of 
motion and which will make the resulting stability analysis tractable. The assumptions 
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are 


(1) The error signals 0 i and 0., i = 1, 2, 3, are available from sensor data and 
can be used for controlling the central body. Letting T^, T c2 , T c2 be the central- 
body control torques about the Xj, yj, and axes of the central body, it is assumed 
that 

T ci = "Vi - C i 0 i i = 1, 2, 3 (4) 

where the and are specified gain constants. 

(2) For controlling the solar array, it is assumed that a sun sensor is used for 
sensing a and that a tachometer or some other rate sensor is used for sensing the 
angular velocity <p of the solar array with respect to the central body. A motor torque 
is used for rotating the solar array relative to the central body. The motor torque T m 
is given by 

T m = K 4 a + C 4^ - a> o ) (5) 

where and are specified gain constants. It is assumed that T m is the torque 
on the central body. 

(3) The satellite remains in a circular synchronous equatorial orbit. 

(4) The angles 0j, 0 2 , 0 3> and 0 4 an< * angular rates 0 0 2 , V an(1 ^4 are 
small. 

(5) The principal axes of inertia of the central body coincide with the central -body 
coordinate axes, and the principal axes of inertia of the solar array coincide with the 
solar -array coordinate axes. The central-body inertias about the x, y, and z axes of 
the central body are denoted by I^ x , Ij y , and I^ z ; and the solar-array inertias about 
the x, y, and z axes of the solar array are denoted by I 2x , I 2y , and I 2z . These in- 
ertias are assumed constant. 

(6) The solar-array center of mass lies along the solar-array shaft. 

(7) The dominant torques on the spacecraft are the control torques given in equations 
(4) and (5) and the constraint torque which constrains the solar array to have only one 
rotational degree of freedom relative to the central body. All other torques are ne- 
glected. 

(8) Nonlinear terms in the equations of motion are neglected. 


EQUATIONS OF MOTION 


The rotational equations of motion are derived in appendix A. The method of equa- 
tion derivation is taken from reference 1. In appendix B, the assumptions of the pre- 
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vious section are used to derive linearized equations of motion. Let u be the vector 
from the central -body center of mass to the solar -array center of mass, and let Uj x , 
u ly, u lz be the components of u in the central-body coordinate system. The linear- 
ized equations of motion can be written in matrix form as 


[M] { 9 } + [C]{ 9 } + [K] { 0 } = { f} 


( 6 ) 


where [M], [C], [K], and {f} are given by 

2 


[M] = 


hx * hx * (I 2z - W 61 " * 
0 


0 


ly 


( I 2 Z " T 2x^ sin ^ cos V 


0 

0 


[c] 


cos <p 

0 

X lz + T 2z ' 

< ! 2z 

- I 2x )sin 2 <p 


0 



X 2y 


0 



V 



‘ 2 2 
u iy + u lz 

” u lx u ly 


" u lx u lz 

0 



M 1 M 2 

_u ly u lx 

u lx + u ly 


-u ly u lz 

0 


(7) 

m t 

- u lz u lx 

" u lz u ly 

2 2 
u lx + u ly 

0 



0 

0 


0 

0_ 



i <f> COS ojj 

0 - w o 

[hx ' T ly + T lz 

+ 2I 2x 

+ %hz- W si " 2 


°1 



-c. 


)[ ! lx * *ly + X lz + 2I 2 x + 2 ' ! 2z ' 0 C 3 " “o[ 2(I 2z " W sin cos <?] 


'4 J 


“ 0^2 


0 

2u ly u lz 

-< 

O' 


2u ly u lz 
„ 2 

0 

2u lx u ly 

0 

(8) 

2u ly 

~ 2u lx u ly 

0 

0 
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[K] = 


K 1 - «o [hz - r lx + X 2x + (I 2z - hJ si » 2 v] 0 

-K^G sin ip 

cos< ?] 

K 4 G sin cp 


'“o[ (I 2z - W sin * cos ?] 

Kg - K 4 -K 4 G cos <p 

0 *3 " “o^lx * *ly + *2x + (I 2z ' ^2x^ cos ^ 


K, 


K 4 G cos cp 



4L-y 

3u lx u ly 

- u lx u lz 

0 

i o,2 Ml M 2 

4u lx u ly 

3 ( u L- u L) 

u ly u lz 

0 

m t 

- 4u lz u lx 

" 3u lz u ly 

2 2 
u lx " u ly 

0 


0 

0 

0 

0 


*4 

0 

K 


4j 


(9) 


{f} = 


m o M 1 M 2 

M™ 


4u ly u lz 

- 3u lx u lz 


u lx u ly 


0 


v. 


( 10 ) 


If the matrices [M], [C], and [K] have constant elements, the stability characteristics of 
the system are found by taking the Laplace transform of matrix equation (6), assuming 
zero initial conditions 


s 2 [M] { 9 } + s[C] { 0 } + [K] { 0 } = { f } (11) 


which can be written as 


[Z(s )]{ 9 } = {f} 


( 12 ) 


where [Z(s)] is given by 


[Z(s)] = s 2 [M] + s[C] + [K] 


We can define a matrix of transfer functions given by 


[Z(s )]" 1 = t A ( s )] 

det[Z(s)] 


(13) 


( 14 ) 
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where 


[A(s)] = adj[Z(s)] 


(15) 


The transformed solution { 6 } is 


{5} = [Z(s)] -1 { f } 


_ r a( s> i { f } 

det[Z(s)] 


(16) 


The solution is obtained by writing the inverse Laplace transformation, 


{ 0 } 


= iT 1 { 0} 


_ jyl/ [A(s)]{ f } \ 
\det[Z(s)] / 


The poles are given by those values of s which satisfy 

det[Z(s)] = 0 


(17) 


(18) 


Equation (18) is called the characteristic equation of the system and the values of s 
satisfying the characteristic equation are called the characteristic roots. If the real 
parts of all characteristic roots are negative, the system is stable. 

If the matrices [M], [C], and [K] are not constant, evaluating the characteristic 
roots by equation (18) is not a strictly valid procedure for determining stability. It is 
customary, however, to evaluate [M], [C], and [K] at a fixed time and to regard the 
characteristic roots of equation (18) as an indication of ’’local” stability at the fixed 
time (ref. 4). 

Precise conditions for which the ’’local" stability procedure is valid are not known. 
It is felt, however, that for slowly time-varying matrices an indication of stability is 
achieved. 

Only those elements of [M], [C], and [K] containing G, sin cp, and cos cp are time 
dependent. The parameter G has a frequency of 1 cycle per year. Terms containing 
sin cp and cos cp have a frequency of 1 cycle per day because cp is nominally equal 


STABILITY FOR u= 0 

We assume in this section that u = 0, or that the centers of mass of the central body 
and solar array coincide. In equations (7), (8), and (9), [M], [C], and [K] are each ex- 
pressed as the sum of two matrices, where one matrix is independent of u and the other 
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is dependent on u. For u = 0, the second matrices on the right-hand side of equations 

(7) , (8), and (9) will be zero. 

Letting be the element in the i^ row and j** 1 column of [Z], equations (7), 

(8) ^ (9), and (13) imply that z^> z j4> z 32> z 34 are zer0 - The characteristic 

equation can be written as 

det[Z] = ( z n z 33 - z l3 z 3i)( z 22 z 44 “ z 24 z 42^ = 0 

The characteristic polynomial contains the term z u z 22 z 33 z 44’ eac ^ z ii> i = 

1, 2, 3, and 4, is a second-degree polynomial in s. Thus, the characteristic equation 
is an eighth-degree polynomial set equal to zero. The characteristic equation is inde- 
pendent of Zgj, Z23, and z^. These four terms are the only ones which contain 
the parameter G. Thus, the stability of the system is independent of the degree of 
coupling of the central-body roll and yaw errors 9 ^ and 9 g into the solar-array point- 
ing error a. 

The factorization of equation (19) implies the separability of the characteristic 
equation into two equations: 


z ll z 33 ~ z 13 z 31 “ 0 ( 2 °) 

z 22 z 44 " z 24 z 42 = 0 ( 21 ) 

Equation (20) contains only roll- and yaw-axis parameters and will be referred to as the 
roll-yaw characteristic equation. Equation (21) contains only pitch-axis parameters and 
will be referred to as the pitch characteristic equation. 


Roll -Yaw Stability 

The roll-yaw characteristic equation is dependent on z^, Zgg, z^g, Zgj. These 
four terms are 

Z 11 = [ (I lx + W + (I 2z " I 2 x )sin V] s2 + [ C 1 + 2w o (I 2z ' l 2x )sin * cos ‘P] 8 

+ K 1 " ^oflz " X 2y + X 2x + ^2z ' W sin2<p ] ( 22 ) 
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z 33 = [ <[ lz + ! 2z> - (I 2z - l2x> sln2 ? , ] s2 + [ C 3 ‘ 2 “o (I 2z ' W 81 " <■> cos *] 8 

+ K 3 - [>lx - ■ly + : 2x + »2z - W cos2 «] < 23 > 

z 13 * [ (I 2z - r 2x )sin V cos 2']= 2 - " 0 [ ! lz - 'lx + 2I 2x + 2(I 2z ' W 8 "' 2 '"] 8 

- “o[< ! 2z - W sin cp cos <pj (24) 

z 3 1 = [ (I 2z - W 8 '" <? cos *>] s 2 - u 0 [li y - I u - l lz - 2I 2x - 2(I 2z - I 2 x )cos 2 ,,] s 

- u 3 [(I 2z - l 2 x' Kin v cos v] (25) 

These four equations can be simplified by using several assumptions. First, assume 
that the inertia difference I 2z - l2 x is much smaller than the inertia sums I lx + I 2x 
and Ij + I 2z . This assumption is satisfied for the satellite configuration shown in fig- 
ure 1. Second, assume that the gain constants Cp C^, Kp and satisfy the follow- 
ing inequalities: 


C i » 2 “o'2x 1 = 1> 3 

Ki » "o ! 2x ‘ - !. 3 

Then equations (22) to (25) can be approximated by 



z U = < I lx + I 2x> 82 +C l 8tK l 


(26) 

z 33 = < ! lz + W® 2 + C 3 8 + K 3 


(27) 

z 13 = 0 


(28) 

Z 31 = 0 

The roll-yaw characteristic equation is then approximated by 


(29) 

[(■lx w® 2 + c i s + K i][ (i iz + y= 2 + c 3 s + K3] = ° 

(30) 
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The first quadratic factor can be thought of as the roll characteristic equation 


( I 1x + I 2x ) s2 + C 1 s + K 1 = 0 (31) 

because only roll-axis parameters appear. Similarly, the second quadratic factor can 
be thought of as the yaw characteristic equation 

(l lz + I 2z )s 2 + C 3 s + K 3 = 0 (32) 

Let s r j and s r 2 be the roots of the roll characteristic equation (31), and let Sy^ and 
Syg be the roots of the yaw characteristic equation (32). Assuming that 

C 1 - 4K i<'ix + W < 0 
C 3 - 4K 3< : lz + ! 2 z> < 0 


then the roll and yaw roots are complex conjugate pairs. 


-C, 


s rl’ s r2 “ 


2 ^lx + W 




^ix + w - c ; 


2 ( I lx + W 


-c. 


s yi’ s y2 - 


^Iz + W 




4K 3«lz + W - C ' 


+ ! 2z> 


(33) 


(34) 


The stability for the roll and yaw axes is determined by the real parts of the roots. 


Re(s rl ) = Re(s r2 ) = 


~ C 1 

^ix + W 


Re(s yl ) = Re(s y2 ) = 


~ C 3 

2 «lz + W 


(35) 


(36) 


Thus, for positive and Cg, the approximate analysis results in stability for the roll 
and yaw axes. If either C-^ or C 3 is zero, there will be roots with zero real parts. 
Thus, proportional feedback alone is not sufficient to achieve stability. 
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A locus of the roots of equations (33) and (34) is given in figure 7. The locus pa- 
rameter is K^y(t)^(I^ x + and the control gains and inertias are specified by 


— = 1 sec i = 1, 3 (37) 

K i 



(38) 


! lx + ! 2 x _ j 
hz + hz 


(39) 


The gain ratios and inertia ratios specified by equations (37) to (39) are chosen to repre- 
sent a realistic dynamic response of a high-power communication satellite configured as 
in figure 1. From equations (37) to (39), the roll roots coincide with the yaw roots in 
figure 7. If equation (37) is replaced by 


C ; 


Kj 


- = 10 sec is 1, 3 


then the roll-yaw root locus is given by figure 8 . With the greater rate feedback, the 
roots are further away from the imaginary axis. 


Pitch Stability 

The pitch characteristic equation is dependent on Zg 2 > 244 , Z 24 , z 42* These 

four terms are 


z 22 = ^y ® 2 + C 2 S + ( K 2 “ K 4^ 

(40) 

z 44 = T 2y s ^ + C 4 S + K 4 

(41) 

z 24 = " C 4 S _ K 4 

(42) 

z 42 = T 2y s2 + K 4 

(43) 
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The pitch characteristic equation (21) can be expanded to yield 




^2y C 2 + I ly C 4 + I 2y C 4^ C 


(Il y K 4 + 


C 2 C 4 + J 2y K 2 


)s* 


+ (C 2 K 4 + C^s + K 2 K 4 = 0 (44) 

Unlike the roll-yaw characteristic equation, the pitch characteristic equation cannot be 
well approximated by two quadratic factors with K 2 and C 2 appearing in one quadratic 
factor and and C 4 appearing in the other. 

One particular case worthy of investigation is the case C 4 = 0. Physically, C 4 
being zero means that no tachometer or other rate sensor is used to measure the body 2 
rate error 9 ^. When C 4 is zero, equation (44) can be expressed as the product of two 
quadratic factors: 

( I 2y s2 + K 4)(V 2tC 2 s + K 2) = 0 (45) 


The characteristic roots dj yK 4 /I 2y , of the first quadratic factor have zero real parts. 
Thus, stability is not achievable for C 4 = 0, and a rate sensor is required to sense 9 ^. 

A locus of the roots of the pitch characteristic equation (44) is given in figure 9. 

The locus parameter is K 4/ /^I 2y , and the control gains and inertias are specified by 


C i 

— =1 sec 
K. 


i = 2, 4 


(46) 


K 4 4 


(47) 


_2y = 1 
r iy " 4 


(48) 


The gain ratios and inertia ratios specified by equations (46) to (48) are chosen to repre- 
sent a realistic dynamic response of a high-power communication satellite configured as 
in figure 1. 

Figure 9 shows that one pair of roots has a positive real part if K 4 /co^I 2y is less 
than about 2xl0 6 . Unlike the roll-yaw case, not all positive values of K 2 , K^, C 2 , and 
C a yield pitch-axis stability. If equation (46) is replaced by 
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'-'i 

— =10 sec i = 2, 4 
K i 

then the locus of the pitch roots is given by figure 10. With the greater rate feedback, 
figure 10 shows that stability is achieved for all positive values of Kg, K^, C 2 , and C^. 
The coefficient of s^ in the pitch characteristic equation (44) is the sum of three 
• terms. If the third term, is dropped, the pitch characteristic equation can be 

expressed as the product of two quadratic factors: 

<‘ly V s4 + < 1 2y c 2 + I ly C 4> s3 + flly K 4 + C 2 C 4 + + <C 2 K 4 + C 4 K 2> S 

+ K 2 K 4 = (l lyS 2 + C 2 S + K2)(l 2y B 2 + C 4 S + K 4 ) (49) 

A coarse approximation to the pitch roots is obtained from the roots of the two quadratic 
factors of equation (49). 


STABILITY ANALYSIS FOR u*0 

Assume now that the central-body and solar-array centers of mass do not coincide 
(u * 0). If { f } in equation (6) is a nonzero matrix, a steady-state error { 0 gg } will re- 
sult. The steady-state error is 


{9 SS 1 = [K]~ 1 { f } (50) 

The matrix { f}, as given by equation (10), will be nonzero only if the vector u has 
more than one nonzero component in the body 1 coordinate system. To evaluate { 0 } 

4 So 

from equation (50), we must first find [K]~ . Assuming that 


» < 


, 4M 1 M 2 1 — I 2 

2z 


i = 1, 2, 3, 4 


the matrix [K] can be approximated as 
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[K] = 


K, 


-K^G sin cp K£ - -K^G cos cp 


K^G sin cp 


K, 


0 

-K, 


K 0 


K^G cos cp 


(51) 


Inverting the [K] matrix of equation (51) and carrying out the matrix multiplication of 
equation (50) yields 


r 


{ ^ss ^ ~ 


M l M 2 “o 
Mr 




4u ly u lz 

K, 


-3ui u, 
lx lz 

K 0 


~ u lx u ly 

K 0 


~ 4u ly u lz G sin ? , 3u lx u lz [ u lx u ly Gcos< ? 


K, 


K r 


K, 


(52) 


The steady-state solar-array pointing error CT gs is calculated from equation (2) to be 


a = 0 
ss 


(53) 


To find the order of magnitude of { 0 }, we assume the following values, which are 

SS 

representative of a high-power communication satellite configured as in figure 1: 


= 450 kg 
M 2 = 150 kg G = 0 


u. = 1 m "' N 


K t = 1 Nm 


cp = 0 


(54) 


The values assumed for u, represent a worst-case situation. For these values { } 

1 SS 

in degrees is 
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{ 9 ssl - 


(55) 


1.4X10" 4 * 

-l.OxlO -4 
-0.3X10" 4 
[l.0xl0“ 4 _ 

Thus, steady-state pointing errors due to nonzero u are negligible. 

To analyze stability for nonzero u, we wish to make assumptions which allow the 
characteristic equation to be separated into three equations for the roll, yaw, and pitch 
axes. Paralleling the stability analysis for zero u, the first step is to factor the char- 
acteristic equation into two quartic factors, where one factor contains roll and yaw pa- 
rameters and the other factor contains pitch parameters. For the zero-u analysis, 
the two quartic factors are given by equations (20) and (21). In order to obtain this fac- 
torization it is necessary for z^ 2 > Zj 4 , Zg 2 > and Zg 4 to be zero. For zero u, this is 
the case as seen by observing the zero elements of [M], [C], and [K] in equations (7), 

(8), and (9). For nonzero u, z^ and z ^ become nonzero. We make the assumption 
that the components of u are small enough to justify setting z^ and z^ to zero. So 
the characteristic equation is still factorable into a roll-yaw characteristic equation, 

z ll z 33 " z 13 z 31 = 0 (56) 

and a pitch characteristic equation. 

z 22 z 44 " z 42 z 24 = 0 

Roll -Yaw Stability 

Continuing to parallel the zero-u analysis, we wish to express the roll-yaw charac- 
teristic equation (56) as the product of two quadratic factors, where one factor contains 
roll parameters and the other factor contains yaw parameters. For the zero-u analysis, 
the factorization is given by equation (30). The required assumptions for making this 
factorization for nonzero u are 


It _ i 

n 1 . 

m,m 9 

1 zz 1 

2x| + 

M t UlxUlz 


« X lx + X 2x 


(58) 
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M 1 M 9 

|*2z " X 2x | + 

m t U1x " 1z 


<<C hz + l 2z 


(59) 


C i » “o 2I 2x + 2 


M 1 M 2 2 


M n 


u iy 


i = 1, 3 


(60) 


K i» < J 2x + 4 


Mj M 2 


u lx u lzj 


i = 1, 3 


(61) 


These assumptions are equivalent to assuming that and z 3 j are zero, and that 

M,Mr 


Z 11 ■ 


T j m l lw 2 (2 2\ 

X lx + r 2x + ~ ~ ( u ly + u lz) 


s' + CjS + K x = 0 


(62) 


z 33 = 


T T M 1 M 2 / 2 2 \ 
I lz + I 2z + -^-( u lx + u lyj 


iy/ 


s‘ + C 3 S + K 3 = 0 


(63) 


Equations (62) and (63) are the separated rqll and yaw characteristic equations, re- 
spectively. 

Define effective inertias 1^ and I z by the equations 


\ = I,„ + Io„ + 


M,M 


l lx + 2x 


12 (2 2 \ 

~mZ W + lz ' 


(64) 


l z = X lz + *2z + 


M 1 M 2 


(-L + u ly) 


(65) 


Assuming that 


C i ' «i^ < 0 


C3 " < 0 


then the roll roots are 
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and the yaw roots are 



( 66 ) 


(67) 


The real parts of the roll and yaw roots are negative for positive and Cg. We con- 
clude that stability is achieved in the roll and yaw axes for nonzero u. Notice, however, 
that the degree of stability is less than for a zero u. By comparing equations (66) 
and (67) with equations (35) and (36), it is seen that the real parts of the roll and yaw 
roots are smaller in magnitude for a nonzero u than for a zero u. For realistic values 
of u, however, the difference in the magnitudes of the roots for the two cases is not 
substantial. The locus of the roll-yaw roots for nonzero u is essentially the same as 
the locus of roll-yaw roots for zero u given in figures 7 and 8. 



then the pitch characteristic equation for a nonzero u differs from the pitch character- 
istic equation for a zero u in that the term 1^ is replaced by an effective pitch inertia 
I'ly defined by 


r ly “ r iy + 


M 1 M 2 


(4 * 4 ) 


From the results of the roll-yaw stability analysis, one might suspect that the ef- 
fective pitch inertia 1'^ will decrease the magnitude of the real parts of the pitch 
characteristic roots and thus decrease the degree of stability. Defining by 
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( 70 ) 





we wish to show that the sum of the real parts of the pitch roots decreases in magnitude 
as AI, increases. If the pitch characteristic equation is normalized so that the coef- 

A y 4 3 

ficient of s is 1, the coefficient of s becomes 


X 2y^ C 2 + + I ly C 4 + AI ly C 4 

( X ly + AX ly) X 2y 


(71) 


Calling this coefficient a, we obtain 


da - l2 2y( C 2 + C 4> 
d < AI ly> (I ly + il ly ) 2 i y 

In terms of the pitch characteristic roots s i = 1 to 4, the coefficient a is 

r 1 

a ' -< s pl + s p2 + s p3 + V 
Since the imaginary parts of the roots sum to zero, 


(72) 


a = -Re(s pl + s p2 + s p3 + s p4 ) (73) 

If it is assumed that C 2 + C 4 is positive, da/d(AI^y) is negative by equation (72). We 
conclude from equations (72) and (73) that the sum of the real parts of the pitch roots 
decreases in magnitude as Al^y increases. We cannot conclude, however, that pitch 
stability decreases as AI 4 y increases. It is possible that the dominant root (the root 
whose real part has the smallest magnitude) is increasing in magnitude at the same time 
that the sum of the real parts of the roots is decreasing in magnitude. Since the degree 
of stability is determined by the dominant root, in such a case the degree of stability in- 
creases as increases. Such a case is shown in the locus of pitch roots shown in 

figure 11. The locus parameter is Al^/I^y and the control gains and inertias are 
assumed to be 

« 2y * r Xy = 400 kg m 2 

4K 4 = K 2 = 40 Nm 
4C4 = C 2 = 40 Nm sec 
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It is seen from figure 11 that the real part of the dominant root increases in mag- 
nitude as Aljy/liy increases from 0 to 0. 7. For realistic configurations of high-power 
communication satellites, Al^/I^y is not greater than 0. 1. The change in magnitude 
of the real part of the dominant root as Al^y/I^ changes from 0 to 0. 1 is not signifi- 
cant. Thus, a nonzero u does not significantly change stability characteristics, and 
the locus of pitch roots for nonzero u is essentially the same as the locus of pitch roots 
for zero u given in figures 9 and 10. 


CONCLUDING REMARKS 

This report presents a small-angle stability analysis for one particular configura- 
tion of a high-power communication satellite having a linear control system. The re- 
sults of the report yield preliminary indications of the relation of stability to satellite 
inertias and control gains. 

Assumptions are made which allow the equations of motion to be linearized and 
which make the resulting stability analysis tractable. The eighth-order characteristic 
equation is separated into two second -order equations for the roll and yaw axes of the 
spacecraft and a fourth-order equation for the pitch axis. For positive control gains, 
stability is achieved for the roll and yaw axes. Not all positive values of control gains 
will yield stability for the pitch axis. Pitch-axis stability is achieved only when the con- 
trol gains are sufficiently large. The degree of stability is not significantly changed 
when the central body and solar array have noncoincident centers of mass. 

The exact roots of the eighth-order characteristic equation were calculated on a 
digital computer by using a subprogram which calculates the roots of a polynomial. The 
exact roots were compared to the roll, yaw, and pitch roots obtained from the approxi- 
mate analysis. For the range of parameters which represent realistic high-power com- 
munication satellites, the approximate roots are within 2 percent of the exact roots. 

Lewis Research Center, 

National Aeronautics and Space Administration, 

Cleveland, Ohio, April 25, 1972, 

164-21. 
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APPENDIX A 


DERIVATION OF GENERALIZED DYNAMIC EQUATIONS 

Introduction 

The generalized nonlinear equations of motion are derived in this appendix. The 
nomenclature and equation derivation is taken from reference 1. The equations as de- 
rived in reference 1 apply to a spacecraft configuration consisting of a central rigid body 
and as many as five auxiliary rigid bodies. Each auxiliary body is assumed to be con- 
nected to the central body at a single point, to have as many as three rotational degrees 
of freedom relative to the central body, and to have no interconnection with other auxil- 
iary bodies. The equations presented in this appendix are a specialization of the equa- 
tions of reference 1 for the case of one auxiliary body (solar array) having one rotational 
degree of freedom with respect to the central body and having the attach point between 
bodies at the center of mass of the auxiliary body. 

Matrix Notation 

It will be convenient to derive equations of motion in matrix notation. Brackets de- 
note a square matrix and braces denote a column matrix. Tor example [A] is a square 
matrix and {u^} is a column matrix. The notation [A] T denotes the transpose of [A]. 

If a. is a vector associated with body i, then { a.j } denotes the column matrix of com- 
ponents of cr expressed in the body j coordinate system. Similarly, { cr. r } denotes 
the column matrix of components of a. expressed in the orbital reference coordinate 
system. Those vectors, such as u, which are not associated with either body 1 or 
body 2 do not have a subscript. For such a vector, { Uj } denotes the column matrix of 
components of u expressed in the body j coordinate system. The components of a 
column matrix are denoted by additional subscripts x, y, and z. Thus, the components 
of the column matrix { } are ay x , a^y, and ct.- z . 

Consider the vector cross product A x B. We wish to express the cross product in 
matrix notation. Let { Aj } and {Bj} be the components of A and B in the body j co- 
ordinate system. Define [J^] to be 

t J Aj l - 


0 

-v 

Ajy 

A. 

0 

-A. 

JZ 

JX 



0 


(Al) 
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Then the components of (A x B) in the body j coordinate system are 

{(AX I).} =[J Aj ]{B j } (A2) 

2 2 

Inertial derivatives of an arbitrary vector a will be denoted by d a/dt, d a/dt ; and 
derivatives relative to a rotating coordinate system will be denoted by a, a. 


Coordinate Systems 

Inertially fixed coordinate system . - The inertially fixed coordinate system will be 
denoted symbolically by 


x 


{Q x } 


Z T 

v. L> 


(A3) 


where Xj, yj, and Zj are unit vectors along the coordinate axes. The origin of the 
{ Qj} coordinate system is at the center of the earth. The vectors Xj and Zj lie in the 
equatorial plane, and Zj points to the sun at the time of the autumnal equinox. The yj 
axis points northward. 

Orbital reference coordinate system . - The orbital reference coordinate system 
will be denoted symbolically by 


r- 

x r 


{ Q r 1 ={y r > 


(A4) 


where x r , y r , and z f are unit vectors along the coordinate axes. The origin of the 
{ Q r } coordinate system is at the center of mass of the central body. The z f axis is 
directed toward the center of the earth, the x r axis is directed along the orbit velocity 
vector, and the y r axis points southward. It is assumed that the satellite orbit is 
exactly the circular synchronous equatorial orbit. 

Body 1 coordinate system . - The body 1 (central body) coordinate system will be 
denoted symbolically by 
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(A5) 


M 


{Qil =<y x > 


l z 


V 


where x^, y^, and z^ are unit vectors along the coordinate axes. The origin of the 
{ Q^} coordinate system is at the center of mass of body 1. The axes of the coordinate 
system are oriented such that they coincide with the corresponding axes of { Q r } when 
body 1 is in its nominal position. 

Body 2 coordinate system . - The body 2 (solar array) coordinate system will be de- 
noted symbolically by 




{Q 2 1 =< y 2 


> 


l z 




(A6) 


where Xg, yg, and Zg are unit vectors along the coordinate axes. The origin of the 
{ Qg} coordinate system is at the attach point connecting body 1 and body 2. The unit 
vector yg is parallel to y ^ and in the same direction as y^ Since body 2 has only 
one rotational degree of freedom with respect to body 1, rotation can take place only 
about the yg-axis. The unit vector Zg is perpendicular to the solar array panels. The 
unit vector Xg is in right-handed orientation to yg and Zg. 


Relative Motion Between Bodies 

Rotational positions . - Transformation matrices from one coordinate system to 
another are [E], [A], and [D] as defined by the following equations. 


{ Q r ) =[E]{Qj} 

(A7) 

{Qj} = [A] { Q r } 

(A8) 

{Q 2 } =[D]{Qj} 

(A9) 


Rotation velocities. - Rotational velocities of coordinate systems are now defined. 
The symbol To is used for inertial rotational velocities, and the symbol S2 is used for 
relative rotational velocities. The following definitions will be used: 
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a> r , fi r inertial rotational velocity of orbital reference coordinate system 

cUj inertial rotational velocity of body 1 coordinate system 

rotational velocity of body 1 coordinate system with respect to orbital reference 
coordinate system 


u >2 inertial rotational velocity of body 2 coordinate system 

fi g rotational velocity of body 2 coordinate system with respect to body 1 coordinate 

system 


The equations of motion will be developed to solve for the time rate of change of the 
relative rotational velocities and fig* 

Rotational accelerations . - Rotational accelerations of coordinate systems are now 
defined. A dot over a column matrix variable indicates time derivative of the column 
matrix components. For example, the three components of { fi rr ) are 



The equations of motion will contain the following rotational accelerations: 

{ d> rr }s{ fi rr ) inertial angular acceleration components of the orbital reference co- 
ordinate system 

{ inertial angular acceleration components of the body 1 coordinate system 

{ } angular acceleration components of the body 1 coordinate system with 

respect to the orbital reference coordinate system 

{ o >22 } inertial angular acceleration components of the body 2 coordinate system 

{ ^ 22 } angular acceleration components of the body 2 coordinate system with 

respect to the body 1 coordinate system 

We derive {cojj} as follows: 

co ^ = fi^, + fi^ (A10) 

{ <X>11 } = + = [A] { ^ rr } + {^ 22 .} (All) 

21> =[ A ](*W + ( ^rr } + {«ii} (A12) 

Letting co Q be the angular velocity of the synchronous orbit, we have 
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(A13) 


'0 ^ 

{ fij. r } = < " W Q 
^0 


{ fi rr } =<0 




Equation (A12) becomes 


We derive { 1 as follows: 


{ <*> 1 1 } = [A] { r } + { ^ n } 


+ S^2 


{ w 22 ^ = { a, 12 ^ + { ^22 ^ “ [®] { } + { ^22 ^ 


(A14) 


(A15) 

(A16) 


Using equation (All), we obtain 


{ ^22 ^ “Plt^H^rr^ + P]{^nl + (^22^ (A17) 

Differentiating equation (A17), we obtain 

W 22 } =[D][A]{S2 rr } +[D][A]{fi rr } +[D]{O n } +[D]{fi n } + {n 22 } (A18) 

In developing equations of motion, we wish to isolate terms containing the relative 
rotational accelerations {fi^} and {0 22 }. We let i u ^ll^ be those terms of {w^} 
which contain { ftjj} and { fi 2 2 } > an< l we let ^n^ll ^ t ^ 10Se terms °f {a> n) which do 
not contain { } and { ft 22 } • We make similar definitions for { u ^ 2 2 ^ and ^n^22 ^ * 

Using equations (A14) and (A18), we have 


{ u u>n} = {ft xl } 

(A19) 

(n^ll^ = 

(A20) 

( u w 22 ^ = [ D ] { ^11 } + (^22^ 

(A21) 

{ n c^ 22 } =[D][A]{fi rr } + [D][A]{n rr } + [D]{fi n } 

(A22) 
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Translational accelerations . - We now seek an expression for the translational ac- 
celeration of the body 2 center of mass expressed in body 1 coordinates. Let u be the 
vector from the body 1 center of mass to the body 2 center of mass. Let S. be the 
vector from the center of the earth to the body i center of mass. Let S m be the vector 
from the center of the earth to the total spacecraft center of mass (fig. 12). We have 
the relation 


M 


1 

u 


S 9 — S + 

* m Mj + M 2 


(A23) 


where and Mg are the masses of body 1 and body 2, respectively. We define G m , 
H u , Gg by the equations 

G m = (A24a) 

dt z 

2 - 

H = — (A24b) 

U At2 



(A25) 


where M^, is the sum of and Mg. Thus, Gg is the translational acceleration of 
the body 2 center of mass. Equation (A25) can be written as 



Mi _ 

— i H 
m t 


u 


(A26) 


The vector H u can be w* 'tten as 


H u = u + 2(a;i x u) + (u> x u) + <x>i x (coj x u) (A27) 

where the derivatives on the right-hand side of the equation are taken with respect to the 
body 1 coordinate system. 

In matrix notation, equation (A27) in the body 1 coordinate system becomes 

<«ul) = i“l) M-Wfu,) +[J lJ11 ]{u l ( +[J u iiP uU ](«i) (A28) 
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Since the vector u has fixed components in the body 1 coordinate system, 


{u^ = {u^ = {0} (A29) 

Equation (A28) simplifies to 

( H ul> = t J iMlK u ll <A30) 

(H ul ) = (A31) 

Separating equation (A31) into two components, { n H u ^ } and { U H U ^ } , 

ln H ul) = -t J ulHn“ll) + [ J ullI J «llH u ll (A32) 

lu H ull =-[ J uJ<u“ll) < A33 > 


Forces and Torques 


In discussing forces and torques, we will use the term "internal torque" to refer to 
reaction torques between the two bodies. All other torques will be termed "external. " 
We make the following definitions: 

F , gravity force on body i, i = 1, 2 

D 

F . external forces, excluding gravity force, on body i 
F r internal reaction force on body 1 due to body 2 

T • gravity gradient torque (about body i center of mass) on body i, i = 1, 2 

& 

T . external torques (about body i center of mass) excluding gravity gradient torque, 
on body i, i = 1, 2 

T r internal reaction torque (about body 2 center of mass) on body 1 due to body 2 
The total force on body 1 in vector and matrix notation is 


F i=F i + F i + F 
el • gl si r 


{ F ell J + {F s11 } + {F r1 } 


gll 


sll 


rl 


(A3 4) 
(A3 5) 


The total torque (about the body 1 center of mass) on body 1 in vector and matrix nota- 
tion is 
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T el = T gl + T sl + < u * F r> + T r I 436 ) 

f T ell I = I T gU ) + ( T sU ) + [J u 1 ][D] t { F r2 } + [D] t ( T r2 } (A37) 

The total force on body 2 in vector and matrix notation is 

F e2 = F g2 + F s 2 ' F r < A38 > 

« F e22> - ( F B 22 > + < F S 22) ‘ < F r2 > < A39 > 


The total torque (about the body 2 center of mass) on body 2 in vector and matrix nota- 
tion is 


T e2 " T g2 + T s 2 ' T r 


(A40) 


^ T e22^ _ i T o-59!^ + " { T v 9 i 


g22 


s22 


r2 


(A41) 


The total force on body 2 is equal to the product of its mass and the acceleration of its 
center of mass 


F e2 = M 2 G 2 


(A42) 


Using equations (A38) and (A42), we have 


F r = F g2 + F s2 ' M 2 G 2 


(A43) 


Body 1 Equation of Motion 

The body 1 equation of motion in vector and tensor notation is 

Tg^ = x Ij • ojj (A44) 

where denotes the derivative of with respect to the body 1 coordinate system, 
and 1^ is the inertia tensor of body 1. 

Expand T e j using equations (A36) and (A43): 

? el = T gl + T sl + u x (F g2 + F s2 - M 2 G 2 ) + T r (A45) 
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(A46) 


T el = T gl + T sl + u X (F g2 - M 2 G 2^ + u X F s 2 + T r 

The term - ^-2^2 is su b rac ti° n of two nearly equal terms. We wish to express 
this term in a more desirable form. Using equation (A23), can be written as 


M l -1 

^ M 2[ S m + TT u l 
V M 2 S 2 _ \ M T 

g2 1-5 13 M. _ 

S + — - u 


(A47) 


Since S m is much larger than (Mj/MrpJu, we can use a perturbation technique to ex- 
press Fgg as 


F g2 = 


-yM, 


3 S 


m 


m 

3 

m 


M i -V 

— u|S. 
M m 


m 


MjU 


m 


M T S 2 
1 m 


(A48) 


We calculate G„, from 
m 


M T G m = £ F gi * E 


r M T S m A - 


i=l 


F . T + 

gi ' Lj si 


L F si 


i=l 


S* i=1 
m 


(A49) 


Equation (A49) reduces to 


g = -I5h + J_ y F . 
m S 3 m t6i 81 

m 


(A50) 


Then MgG 2 is obtained by using equation (A26) 


_ _ M.M 9 _ 

M 9 G 9 = M 9 G rn + — — - B 
2 2 2 m u 


(A51) 
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Now using equation (A50), 


m 2 g 2 = " 


^ M 2 S m , M 2 f F , M 1 M 2 „ 

0 3 M T S1 Mrp U 


(A52) 


m 


Subtracting equation (A52) from equation (A48), we obtain the term - MgGg) in a 
more desirable form. 


F g2 - M 2 G 2 = ' yM 2 


“i' 5 .-*. , M i“ 


M T S m 


m t4i 


M l M 2i 

m T^i S1 M t ^ 


(A53) 


Combining equations (A46) and (A53) yields 


T el - T gl + T sl + u x 


3xM l M 2< S m ' M 2 y ? M l M 2 H u 

M T S m Mt 1=1 S1 Mt 


+ u x F s2 h- T r 


(A54) 


Using the vector identity 


<“• V<" x V = “ x [ S m x < S m x “>] 


(A55) 


equation (A54) becomes 


M 2 _ 

T- 1 ~ T~. + T„. + u X lF g 2 - — — ^ F si |+T r 

m T i=l 


el gl + si 


3y Ml M 2 {u X [s m X (S m X u)jj 


M t M 2 (u X H u ) 


S m M T 


M n 


(A56) 


We now wish to write the vector equation (A56) in matrix notation. The components of 
the vector S m in the orbital reference coordinate system are given by 
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0 


< S mrH 0 


where is the magnitude of S . Then 
m m 


< S ml> =- S mt A > 0 - S m ^ A 3> 


where { A^} denotes the third column of the matrix [A]. We can write equation (A56) in 
the body 1 coordinate system as 

t T ell) - (TgU 1 + 1T sll> + t J ul](~ tF Sll » + ^ < F s2ll) 

3yM 1 M 0 M.M, / \ 

- -r- 5 -* [JulI J A3][ J A 3 H “1 1 - [ J ul) (i n H ul) + ( u H ul l) (*») 


4 M T 


Rewriting equation (A44) in matrix notation. 


f I l^( { n cb ll } + ^11*) + [ J collW^ll} " {T eir 


(A60) 


where [1^] is the inertia matrix of body 1 with respect to the body 1 coordinate system. 
Combining equations (A59) and (A60) 


“i “ 2 


PlHu“ll> + ^l J ulllu H ull = (In l + (T rl ) 


where {q^} is defined to be 
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{ ^11 } - 'IWn^ll* " + {T gll } + {T sll* 

/ M 9 M« \ 3yM 1 M 0 

+ [ J ul]("“ {F sll } + ~ ^ F s 21 } ) + ~T [ J ultf J A3tf J A3] {u l } 

\ M T M T / SlM™ 


M t M 0 

( A62 > 

By using equations (A19) and (A32), equation (A61) becomes 

[y< n U> - — 5 [ J ul]f J ulH n ll } - did * ( T rll < A63) 

In the analysis to follow, it will be convenient to express equation (A63) as 

[ L ll]{^ll} + [^"T 2 ]{^ 22 } = + ^rl^ (A64) 

where [L 1;L ] and [L 12 ] are given by 

M,M 9 

t L ll] = Ull " [ J ulH J ul] (A65) 

[L 12 ] = [0] (A66) 


Body 2 Equation of Motion 

The body 2 equation of motion in vector notation is 

T g 2 — I 2 * a? 2 + u> 2 ^ I 2 * c^2 (A67) 

where ci 3 2 denotes the derivative of w 2 with respect to the body 2 coordinate system 
and I 2 is the inertia tensor of body 2. From equation (A40), T g2 is given by 

f e2 = T s2 + T s2 - T r (A68) 

Combining equations (A67) and (A68) in a matrix equation in the body 2 coordinate system, 
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P2l(ln w 22J + ^ u w 22 + t J cu22^ I 2^ w 22 ^ ~ ^ T g22^ + {T s22^ " {T r2* < A69 ) 

where [I 2 ] is the inertia matrix of body 2 with respect to the body 2 coordinate system. 


By using equation (A21), equation (A69) becomes 

+ {^22^) = ^22^ ~ ^r2^ (A70) 

where { q 22 } is given by 

^22^ = n^22^ " ^0)22^2^^22 E + ^ T g22^ + ^ T s22^ ( A71 ) 

Equation (A70) is rewritten as 

[L 2 i]{fii 1 } +[L 22 ]{fi 2 2^ = ^ q 22^ ” ^ T r2^ (A72) 

where [L 2 j] and [L 22 ] are given by 

[L 2 i3 = [I 2 ][D] (A73) 

[L 22 ] = [I 2 ] (A74) 


Elimination of Locked-Gimbal Reaction Torques from Equations of Motion 

We repeat the body 1 and body 2 equations of motion 

+ [Lj 2 ]{ ^ 22 } = {q^j} + [ D ]'i'{ T r 2} (A75) 

[ L 2l]{^nl + t L 22^^22^ = ^ q 22^ " ^ T r2^ ( A76 ) 

Body 2 has rotational freedom only about its y-axis. The x-axis and z-axis of body 2 
are locked gimbals. Then t T r2> is composed of two parts 

( T r2 1 = ( T fg2 1 + i T lq2 } (A77) 

where { T F q 2 } is the reaction torque about the "free-gimbal" y-axis, and { T L q 2 ) is 
the reaction torque about the '’locked-gimbal" x-axis and z-axis. {T FG2 } has only a 
y component, and { T^^ } has only x and z components. 
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0 


^ T FG2 ^ - 


^ T LG2 ^ ” 


FG2y 

L 0 . 

LG2x 

0 


L T LG2z 


(A78) 


(A79) 


Substituting equation (A77) into equations (A75) and (A76), 

[Ln]{fin} + [L { ^22 ^ = ^11^ + [^1 t^FG2^ + ^]t^LG2^ (A80) 

+[L 2 2]{^22^ = ^22^ “ ^ T FG2 ^ " ^ T LG2^ (A81) 


We now manipulate the body 1 and body 2 equations (A80) and (A81) in order to 
eliminate { T LG2 } . It will be convenient to introduce three new matrices, [U], [P], 
and [U - P], defined by 


[U] = 


[P] = 


[U - P] = 


0 

1 

0 

0 

1 


0 0 

1 
0 
0 


(A82) 


(A83) 


0 

0 

0 


(A84) 


We have the following relation: 
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FG2 

LG2 


} " ^ T FG2 ^ 

[U- P]|T FG2 ) = { 0} 

(A85) 

} = {0} 

[U - P]{ t lG2^ * ( T LG2> 

(A86) 


To eliminate { t lq 2 ^ ^ rom the body 2 equation of motion, we premultiply equation (A81) 

by [P]. 


[ p ][ L 2iK ^11 } + [ p ][ L 22^ ^22 } _ [ p ]i c l22^ ~ ^ T FG2^ 


(A87) 


To eliminate { T^q 2 } f rom the body 1 equation of motion, we must first premultiply the 
body 2 equation (A81) by [D] T [U - P]. 

[D] T [U-P][L 21 ]{fin} +[D] T [U-P][L 22 ]{fi 22 } =[D] T [U-P]{q 22 } 

(A88) 


By adding equation (A88) to the body 1 equation (A80), we obtain the new body 1 equation 
with { T L q 2 } eliminated. 

([Lll] + P] T [0-P][L 21 ]){O U > + ([L 12 ]*tD] x [U-P][L 22 ])(Si 22 ) = {4 U ( 


• [D] t 1 T FG2 + [D X ][U - P] ( q 22 ) (A89) 

Rewriting the new body 1 and body 2 equations (A89) and (A87), 

[M n ]{fi u } +[M 12 ]{6 22 } = {b u } (A90) 

t^21^ ^11 ^ + [•^'22^ ^22 ^ = { ^22 ^ (A91) 

where [M u ], [M 12 ], [M 21 ], [M 22 ], {bill, and {b 22 } are given by 

[M n ] = [L n ] + [D] t [U - P][L 21 ] (A92) 

[M 12 ] = [L 12 ] + [D] t [U - P][L 22 ] (A93) 

[M 21 ] = [P][L 21 ] (A94) 

[M 22 ] = [P][L 22 ] (A95) 
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{b u } ={q u } +[D] t {T fg2 } +[D] T [U-P]{q 22 } 


(A96a) 


{b 22 } =[P]{q 22 } - { T FG2 } (A96b) 

Since [M^], [M 22 ], {b 22 } all involve a premultiplication by [P], the first and third 
rows of these three matrices are zeros. Also, the first and third components of { fi 22 ) 
are zero because body 2 has rotation only about its y-axis. Combining equations (A90) 
and (A91) into a single matrix equation, we have 



(A97) 


The third and fifth rows and columns of the 6X6 coefficient matrix can be deleted to form 
a 4X4 system 

M 12xy ^llx 
M 12yy ^lly 

< • r 

M 12zy fi llz 
M 22yy. [ fi 22y^ 

For shorthand notation, equation (A98) will be written as 

= {b } (A99) 


' M llxx 

M llxy 

M llxz 

M llyx 

M nyy 

M llyz 

M llzx 

M llzy 

M llzz 

. M 21yx 

M 2lyy 

M 21yz 
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APPENDIX B 


LINEARIZATION OF DYNAMIC EQUATIONS 
Assumptions 


The dynamic equations derived in appendix A will now be linearized by assuming 
that the angular coordinates 9 ^ are small. Thus, sin 0^ is approximated by 9 ^, and 
cos 9 . is approximated by 1. We make the following assumptions: 

(1) The body 1 coordinate axes are nearly the same as the orbital reference coor- 
dinate axes. Thus, we can specify the orientation of the body 1 coordinate system with 
respect to the orbital reference coordinate system by three small angles 0j, 9 2> 
and 0g. The orbital reference coordinate axes, if rotated through an angle 9 ^ about 
the x r -axis, through an angle 0g about the y r -axis, and through an angle 9 3 about 
the z r -axis, will coincide with the body 1 coordinate axes. The transformation ma- 
trix [A] becomes 



"3 

1 




(Bl) 


(2) The body 1 and body 2 inertia matrices are diagonal matrices. Equivalently, the 
geometric axes coincide with the principal axes of bodies 1 and 2. The inertia matrices 
[IjJ and [Ig] become 


I 


[Il] = 


lx 

0 


I 


Pal- 


2x 

0 

0 


I 


iy 

0 


1, 


2y 

0 


0 

0 

^z. 

0 

0 

r 2z 


(B2) 


(B3) 


(3) The body 1 and body 2 inertia matrices have constant elements, and the body 1 
and body 2 masses are constant. 
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(4) The dominant torques on the spacecraft are the control torques. Environmental 
torques such as the gravity gradient, solar pressure, and magnetic torques are negli- 
gible. Any friction torque between the two bodies is negligible. It is assumed that 0- 
and i = 1, 2, 3 are known from sensor data and that the components of the body 1 
control torque vector T c consist of proportional and derivative feedback. The vector 
T c is given in the body 1 coordinate system by 


r - K ih - c iV 


l T cl> * V K 2 s 2 ■ C 2 e 2 > 


^ -K 3 0 3 " C 3 0 3-> 1 


(B4) 


(5) A motor torque is used for rotating body 2 with respect to body 1. See figure 3 
for the definition of the angle cp . The angular velocity of the body 2 coordinate system 
with respect to the body 1 coordinate system is cp. The nominal value of cp is one 
earth rate w 0 . So the rate error is given by cp - u) Q . It is assumed that the motor 
torque responds to the position error a, as defined in figure 4, and the rate error 
cp - The motor torque vector T m is given in the body 2 coordinate system by 


r 


o 


*\ 


i T m2>=l¥ + C #- w o)[ 
0 


v. 


(B5) 


where it is assumed T_ acts on body 1. 

m J 

(6) Nonlinear terms in the equations of motion are negligible. Thus, terms such as 

2 

9 . , and 9^9^ are removed from the equations of motion. 


Calculation of o 

The fourth angular coordinate 9 ^ is defined as 

£4 = <P ~ (a + a > Q t) (B6) 

where the angles a and a) Q t are shown in figure 5. 

The body 2 pointing error a is equal to 9 ^ when 9^ = 9^ = $3 = 0. When body 1 
pointing errors are nonzero, they contribute to o. Let 6 be a unit vector directed from 
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the center of the earth to the sun. Let 6p be a vector formed from the projection of 6 
into the Zj-Xj plane. Let 0 g be the angle between Zj and 6, and let 0 p be the angle 

between z T and 6 . The geometry is shown in figure 13. 

A r 

The angles 0 O and 0 are related by 
s p 

tan 0_ cos i = tan 9 (B7a) 

s p 

where i is the inclination of the ecliptic plane to the equatorial plane (i = 23. 449°). In 
the inertial coordinate system, 6 is given by 


''sin 0 cos O 
s 


tty H 


-sin 0 sin if 
s 

cos 0„ 


(B7b) 


A series of transformation matrices is required to express 6 in the body 2 coordinate 
system. 


{6 2 } = [D][A][E] { 6j] 


The [A] matrix is given by equation (Bl), and [D] 


[D] = 


’cos (p 
0 

J_sin cp 


[E] = 


-cos(0 + a 2 + u> 0 t) 


sin(0 + a 1 + co Q t) 


0 

1 

0 

0 

-1 

0 


and [E] are given by 

-sin <p" 

0 

cos cp m 

sin(0 p + otj + u> Q t) 
0 

COS(0 p + QTj + C0 o t) 


(B8) 


(B9) 


(BIO) 


The angle a is assumed small, and the sign convention for a is that it is negative if 
the projection of 6 into the ^2~ x 2 P lane lies between the z 2 ~axis and x 2 -axis. Thus, 


a = 


sin a = 



(Bll) 
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Using equations (B 6 ) to (Bll), cr is calculated to be 

o = Jg sin(<p - ^ 4 )]^! + 0 2 + g|cos(^ - ^ 4 )]^3 + ^4 
where G is given by 


-sin 0 g sin i 
^1 - sin^Sg sin^i 


(B12) 


(B13) 


Since 9 4 is a small angle and cp in general is not small, we can approximate equa- 
tion (B 12 ) by 


a = (G sin cp)Q ^ + (G cos cp)9 g + 9 4 

The motor torque of equation (B5) can now be written as 


t T m 2 } H (K 4 G sin <p) 0 j + K 4 9 2 + (K 4 G cos <p)9 3 + K 4 0 4 + 




(B14a) 


(B 14b) 


Calculation of CM3 

The matrix [M] of equation (A99) consists of elements of the matrices [M^], [Mj 2 ], 
[M 21 ], and [M 22 ], which in turn consist of elements of the matrices [Ljj], [Lj 2 ], [Lgj], 
and[L 22 ]. Evaluating [L^], [L 12 ], [L 21 ], and [L 22 ] from equations (A6 5), (A 66 ), (A73), 
and (A74), 


[L U ] = 


iix 

0 

0 



2 

u lz 

~ u lx u ly 

' u lx u lz 

0 


0 

M 1 M 2 


2 2 



+ 

" u ly u lx 

u lz + u ly 

" u ly u lz 
2 2 

0 

0 

*lz_ 


_ - u lz u lx 

~ u lz u ly 

u lx + u ly 


(B15) 




'0 

0 

Lo 


0 0 

0 0 

0 0. 


(Bie) 
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t L 2 l] = 


I 2x cos cp 


I 2z sin cp 


I, 


2y 

o 


[ L 22-l = 


2x 

0 


2y 

0 


_I 2x sin 
0 

l 2z cos cp 

0 
0 


I 


2z 


Evaluating [M^ ^], [M 12 ], [M 21 ], and [ M 22^ from equations (A92) to (A95), 

.2 


[“ni - 


! lx + : 2x + ( hz - W sln <" 

0 I 

«2z- l 2 x )sin cp cos cp 


iy 


( ! 2 Z - I 2x )sin cp cos cp 
0 

; _2 


0 l lz + l 2z - ( X 2z - W sin 


M 1 M 2 


M n 


2 2 

' u lx u ly 


ly + u lz 

-u lx u lz 

" u ly u lx 

2 2 
u lz + u ly 

“ u ly u lz 

,- u lz u lx 

_u lz u ly 

u lx + u ly 


[M i2 ] = 


2x cos 

0 

J 2z sin V 


0 

0 

0 

(B20) 

X 2x sin 

0 

X 2z cos 


ro 

0 

ol 



t M 2l] = 


I, 


2y 

o 


0 

0 . 


(B17) 


(B18) 


(B19) 


(B21) 


[ m 22 ] = 


0 

LO 


I, 


2y 

o 


o 

a 


(B22) 
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When [M^], [M^]* [ M 21^ and [ M 22^ are combined into a 6*6 matrix (see eq. (A97)) 
and the third and fifth rows and columns are deleted, the matrix [M] is obtained 


X lx + X 2x + ^2z " W sm ^ 


[M] = 


0 

(I 2z ' I 2x ) sin <P cos <p 
0 



0 



«2z- l 2 x )sin <p cos (p 
0 

T lz + X 2z “ ^2z ” W sin2< P 
0 


0 

0 

0 

l 2y. 



' 2 2 
U lx + u ly 

” u lx u ly 

' u lx u lz 

0 

+ M 1 M 2 

“ u ly u lx 

2 2 
u lz + u lx 

~ u ly u lz 

0 

m t 

" u lz u lx 

" u lz u ly 

u lx + u ly 

0 


0 

0 

0 

0 


(B23) 


Calculation of {Q} 


The first three components of the column vector { & } of equation (A99) are the com- 
ponents of the vector in the body 1 coordinate system. Because of our linearity 
assumptions, these components are simply Q^, ^ and The fourth component of 
{ fi} is the angular acceleration <p of body 2 with respect to body 1. From equa- 
tion (B6), we can solve for cp in terms of 9^. 

<p = (B24) 


Thus, { n } reduces to 




9 

w =< .. y 

e 




(B25) 
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Calculation of {b} 


In calculating the {b} matrix, we neglect all nonlinear terms. The {b} matrix of 
equation (A99) consists of components of the matrices {b^} and { bg 2 } • As a 
step to calculating { b^ } and { b£ 2 ) ^ rom equations (A96) and (A97), we calculate { q^ } 
and { q 2 2 } from equations (A62) and (A71). Our assumption that forces and torques 
other than control forces and torques are negligible implies that 


{ T gii > = { F sii> ={°> i = 1 > 2 


So {q 11 } and { q 2 2 } reduce to 


M 1 M 2 






3yM,M 9 

+ < T sH> + “ [ J ul]l J A3^ J A3^ u l } 


S m M T 


(B26) 


^22^ ~ 't J 2^ n w 22^ " t J oo 22^*2^ w 22 ^ + ^ T s22^ 


(B27) 


We now proceed to evaluate each of the five terms on the right-hand side of equa- 
tion (B26). 

The matrix { n co ^ i } , from equations (A20), (Bl), and (A13), is 



p 


■ 


r 


r 


0 

h 



0 


_a, o 0 3 

=lA]fn rr ) - 

"®3 

0 

h 

< 

00 

0 


0 ► 



-®1 

0 


0 

^ J 


00 

v. o 1 J 


The first term of equation (B26) is then 


fli „co 0, ^ 


-IWn^ll} H 


lx o 3 
0 


L^lz^o 0 !, 


(B28) 


(B29) 
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The matrix {co^}, from equation (All), is 


{ ^ 1 1 } = [A]{ fi rr } + { ^ii } - 


rr J 


‘11 


1 

-dr. 


- 9 , 


r o^ 


< w < 

0 
v j 


V 


r + V2 

V.V 


r V u o 9 3 


{ } H 0 5 " W 


The second term of equation (B26) is then 


^3 +w o 01 


-[<WW w u> =1 


"lly^llz^lz 

^llz^llx^lx 


- V] K^3 + w o 0 l^ r lz - hy 
: w o ( *l + w o*3 )(I ly " X lx 


C’llX^lly^ly “ hx>J 


The matrix { n H ul }, from equation (A32), (B28), and (B29), is 
{ n H ul } = ~[ J ul^n w ll } + [ J «llF«llH u l> 


in H ul} H 


- 2w 0 u 2*l + 2w o u 1^2 + w o u 2 0 3 " u l w 5 j 

_ A . u. r 


-Ci) Uo 0 1 + O) Ui0 o 

o 3 1 o 1 3 


v 2w o u 3®2 - 2w o u 2®3 " w oVl " u 3 w o J 


The third term of equation (B26) is then 


(B30) 


(B31) 


► (B32) 

i 

(B33) 

(B34) 
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MjMg 

Mrri 


r 2a) o u 3 u 2^2 - 2 a, o u 2^3 ^ 


Mi m 2 


[ J ulKn H ul}=~ 


Mr 


- 2w o u 2 u 3 9 1 + 2 a, o u 2 u l 0 3 
2 ct) o u Pl “ 2 w o u 1 u 2^2 


r 9 ? 90 

cAofl, - O)„Uo0 


M 1 M 2 


o 3 1 


o u 2 u l " w o u l u 3 0 3 


w o u 1 u 2 9 1 + w o u 2 u 3 S 3 


2 2 2 2 2 

"CO U-«U<j5 1 - (jjAlcydn + 0) U,0o 

0131 o 2 3 ol3 


(B35) 


The fourth term of equation (B26) consists of the control torque of equation (B4). 


i T s ii> -< 


-KjOi - C 1 S 1 ' 
“ K 2 0 2 ‘ C 2®2 


v. K 3 e 3 " C 3 0 3y 


(B36) 


The matrix { Ag} is the third column of [A]. 


r-e: 


{Ao} =< 




(B37) 


The fifth term of equation (B26) is then 


SyMiMq 3yMrM 9 


S m M T 

m i 


S m M T 


u 2 9 1 " u 3 0 1 ‘ U 1 U 2 0 2 + U 2 U 3^ 
<~ u 1 u 2 9 1 + u 1 9 2 ” u 3 0 2 “ u l u 3 


U 3 U 1 0 1 + U 3 U 2 0 2 


(B38) 


The matrix {q^} is then given by the sum of equations (B29), (B32), (B35), (B36), 

and (B38). 

To evaluate { q 2 2 ^ we eva ^ ua ^ e each of the three terms on the right-hand side of 
equation (B 27). The matrix { n ^ 22 }, from equations (A22), (B9), (Bl), and (A13), is 
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r 


o 


{ n d> 22 ^ = ,< ^co 2 cos cpj 9 j + (2co o cos cp)d^ - (w 2 sin cpj 9^ + (2w Q sin cp)9 ^ 
-(co 2 sin c pjd^ - (2o? 0 sin < p)6^ - cos <p)<9g + (2cu 0 cos cp)9^J 


(B39) 


The first term of equation (B27) is then 


-ty { 11^22 ^ = * X 2y [" ( w o cos <?) 5 1 " ( 2a> 0 cos ^3 + ( w o sin *?) 9 3 " ( 2a, o sin ^l] 

J 2z [( w o sin <?) 9 1 + ( 2a, 0 sin ^)®3 + ( w o cos <?) 0 3 " ( 2co 0 cos ^)®l] J 


(B40) 


Using equations (A16), (B31), (B9), and (B3), the second term of equation (B27) is found 
to consist entirely of nonlinear terms. 


-lWWi w 22> = <°> ( B41 > 

The third term of equation (B27) is also { 0 } . 


{ T S 22 } = { ( B42 ) 

The matrix { q22 ^ is then given by the sum of equations (B40), (B41), and (B42). 

Knowing {q^} and { ^22 ^ we can now calculate {b^} and { ^ 22 } * rom equations 
(A96a) and (A96b). Recalling that { Tp^ ) is the motor torque of equation (B14b), 


r 


- ^ 


{ t fG 2 } =< { ( K 4 G sin <?) 0 i + K 4 0 2 + ( K 4 G cos ^ d 3 + K 4 0 4 + C 4^4 

0 J 


we can calculate {b} in the matrix form 

{b} = -[C]{ 0} - [K] { 9 } + {f } (B43) 

where [C], [K], and {f} are given by 
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[C] 


C 1 + “o 2(I 2z ' l 2 x )sin cp cos 0 0 

0 c 2 

*0 fix * T ly + X lz + 2I 2x + *hz ‘ W 008 ^] 0 

0 0 


'"o fix - J ly + ! lz + 2I 2x + 2 »2z ' W 81 " 2 J 0 


C 3 * “ 0 2(I 2z ' W 8 ™ V cos cp 
0 


“o M 1 M 2 


Mp 


" 2u ly u lz 


2u 


iy 


2u ly u lz 


" 2u lx u ly 


-2u 


iy 

2u lx u ly 

0 

0 


- C 4 

0 

0 

0 

0 

0 

0 


[K] = 


k 2 -k 4 


^-“otlz^lx + ^x^^z-W 8 ^] 0 

-KjG sin cp 

- w o[ (I 2z - W sin< ? cos< ?] 

K 4 G sin cp 


- w o[ (I 2z ' r 2x )si " V cos <p] 

-K 4 G cos cp 

K 3 ' w o [ J lx ' r ly + r 2x + (I 2z ' W 008 ^] 
K 4 G cos cp 


co 2 M,M, 
o 1 2 




l ( u lz' u ly) 


4u lx u ly 

- 4u lz u lx 


-K„ 


3u lx u ly 

’Mz'O 


- 3u lz u ly 


- u lx u lz 0 

u ly u lz 0 

Ax - Ay 01 


r4u ly u lz 


{f} = 


w o M 1 M 2 


M n 


< 


" 3u lx u lz 


u lx 11 ly 


0 


v. J 

By using equation (B43) for {b} , the final linearized equations of motion can now be 
written from equation (A99). 


[M ]{ 9 } + [C]{0} + [K] { 0 } = { f} 


E-6865 
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APPENDIX C 


SYMBOLS 

coefficient of s in pitch characteristic equation 
coefficient matrix of { 6 } 

i in control gain for rate feedback, i = 1, 2, 3, 4 
torque matrix 

parameter determining degree of participation of and 9 g in ex- 
pression for cr 

effective satellite roll inertia 

effective satellite yaw inertia 

mass moments of inertia of central body about x, y, z axes of 
central-body coordinate system 

mass moments of inertia of solar array about x, y, z axes of solar- 
array coordinate system 

effective central-body pitch inertia 

variation of I' ly from I^ y 

inclination of equatorial plane to ecliptic plane 

coefficient matrix of { 9 } 

i** 1 control gain for position feedback, i = 1, 2, 3, 4 

line formed by projecting earth- sun line into equatorial plane 

coefficient matrix of { 9 } 

sum of masses of central body and solar array 

mass of central body 

mass of solar array 

real part of 

Laplace operator 

iu 

i in root of pitch characteristic equation, i = 1, 2, 3, 4 

i th root of roll characteristic equation, i = 1, 2 
til 

i tn root of yaw characteristic equation, i = 1, 2 



T cP T c2’ T c3 


m 


u 


u lx’ u ly> u lz 
x p y p ij 

V V i r 

V V K 

x i> yp z i 

Xp yp z\ 

x 2’ ^2’ z 2 
[Z] 

z ij 


y 

6 

*i 

0 s 

{0} 

{0} 

{$} 

^ 0 ss^ 
a 

<P 


central-body control torques about x, y, z axes of central-body co- 
ordinate system 

motor torque on central body 

vector from central-body center of mass to solar-array center of 
mass 

x, y, z components of vector u in central-body coordinate system 

unit vectors along axes of inertially fixed coordinate system 

unit vectors along axes of orbital reference coordinate system 

unit vectors along axes of orbital reference coordinate system with 
origin translated to center of earth 

unit vectors along axes of central-body coordinate system 

unit vectors along axes of central -body coordinate system with origin 
translated to solar-array center of mass 

unit vectors along axes of solar-array coordinate system 

impedance matrix 

1U i,L 

element in i row, j n column of [Z] 

angle subtended by line L and vector z' r evaluated at time equal to 
zero 

earth gravitational constant 

unit vector directed from solar-array center of mass to sun 
unit vector formed by projecting 6 into equatorial plane 

it. 

i angular coordinate, i = 1, 2, 3, 4 

apparent angular displacement of sun from autumnal equinox measured 
in ecliptic plane 

column matrix of angular coordinates 
column matrix of angular rates 
column matrix of angular accelerations 
column matrix of steady- state displacements 
solar -array error 

angular rotation of solar-array coordinate system with respect to 
central-body coordinate system 
angular velocity of circular synchronous orbit 
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Figure 1. - Configuration of high-power communication satellite. 
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Figure 6. - Sun angle 9 S measured in ecliptic plane. 


54 
















NATIONAL AERONAUTICS AND SPACE ADM ISTRATION 
WASHINGTON. D.C. 20546 


OFFICIAL BUSINESS 
PENALTY FOR PRIVATE USE $300 


FIRST CLASS MAIL 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 



NASA 451 






POSTMASTER: 


If Undeliverable ( Section 158 
Postal Manual) Do Not Return 


'The aeronautical and space activities of the United States shall be 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof.” 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. 

CONTRACTOR REPORTS. Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include conference proceedings, 
monographs, data compilations, handbooks, 
sourcebooks, and special bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 


